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LONG-TERM  BIOEFFECTS  OF  435-WZ  RADIOFREQLENCT  RADIATION 
ON  SELECTED  BLOCO-BOWC  ENDPOINTS  IN  CANNULATED  RATS 
VoIum  3.  Plasma  Prolactin 

I.  INTRODUCTION 

Throughout  the  developed  world*  and  particularly  in  the  United  States*  the 
20th  century  has  marked  a  period  of  tremendous  progress  In  communication, 
information,  and  electronic  sciences.  Many  of  the  major  technological  advances 
during  this  period  Involved  transmitting  energy  over  vast  distances  using 
electromagnetic  waves.  This  progress  had  the  side  effect  of  altering  the 
planet’s  electromagnetic  environment.  Radio*  radar*  and  television 
transmissions  have  increased  the  ambient  electromagnetic  radiation  level  by 
several  orders  of  magnitude.  At  this  time*  despite  many  studies  performed  in 
this  field,  the  biological  effects  of  this  omnipresent  el ectromagnet ic 
environment  on  organisms  are  not  well  understood. 

This  report  presents  results  of  plasma  prolactin  levels  measured  in  200 
male  Sprague-Oawley  rats  chronically  exposed  to  a  1.0  mW/crir^,  435-MHz  pulsed- 
wave  (1.0  us  pulse  width,  1-kHz  pulse  rate)  electromagnetic  environment  for  a  5- 
month  duration.  The  exposure  group  consisted  of  ICO  cannulated  rats  housec  Jn 
Plexiglas  cages  arrayed  on  the  tiers  of  a  stacked*  parallel-plate  circular 
waveguide.  Engineering  aspects  of  this  waveguide  and  the  exposure  environment 
it  generated  have  been  previously  reported  [1].  The  sham-exposure  group 
consisted  of  100  cannulated  rats  housed  in  an  identical*  out  unenergized, 
collocated  facility.  The  biological  effects  of  this  radiof requency  radiation 
(RFR)  exposure  on  plasma  adrenocorticotrop 1c  hormone  (ACTH)  and  plasma 
corticosterone  concentrations  in  the  same  animals  have  already  been  recorded 
[21. 

3rolactin  was  identified  50  years  ago  as  a  lactogenic  hormone  secreted  oy 
the  anterior  pituitary.  Recently  developed  sensitive  and  specific 
radioimmunoassay  methods  have  led  to  knowledge  of  the  physiology  and 
oatnophysiology  of  prolactin  secretion.  3rolactin  is  apparently  releasee  -'n  a 
pulsatile  fashion  [3].  The  pulses  are  small*  except  during  sleep  when  rnarxec 
rises  In  prolactin  concentrations  have  been  noted.  Plasma  prolactin  level  in 
undisturbed  intact  male  rats  Is  about  10-15  ng/mL  '•*].  In  male  animals,  p'asma 
prolactin  levels  rapidly  Increase  from  4  to  10  times  the  oasal  level  in  response 
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to  various  stressors  [5,6].  Some  known  stressors  Include  surgery  or  anesthesia 
(Increase  or  decrease  [7]),  feeding  [8],  and  brief  handling  or  mild  ether 
exposure  [4,9].  Elevated  ambient  temperature  (36  °C  for  20  to  360  min) 
associated  with  body  hyperthermia  evokes  Increases  In  circulating  levels  of 
prolactin  CIO],  During  stress,  prolactin  Is  released  In  a  quantitative  fashion 
[11];  thus,  the  level  of  plasma  prolactin  can  be  used  to  measure  the  level  of 
stress  [12]. 

Intraventricular  brain  Injection  of  E-endorphin  (£  L?Hgi_g^)  1n  urethane- 
anesthetized  male  rats  leads  to  a  dose-dependent  Increase  of  plasma  prolactin 
levels  [13].  Thus,  the  plasma  prolactin  level  corresponds  to  the  level  of 
stress  In  a  fashion  similar  to  plasma  ACTH  and  plasma  corticosterone.  Male  rats 
acutely  exposed  to  visual  or  audiogenic  stimulation  exhibit  rapid  and  marked 
prolactin  secretory  responses  [14].  This  suggests  that  the  response  to  an  acute 
exposure  of  neurogenic  stress  In  the  male  rat  Is  elicited  via  a  neural  pathway 
Impinging  upon  the  medial  basal  hypothalamus  from  the  rostral  direction. 
Midbrain  lesions  slightly  alter  the  level  of  plasma  prolactin  In  adult  male  rats 
[IS],  but  the  Integrity  of  the  amygdala  Is  not  essential  for  the  normal  basal 
and  diurnal  hormone  profile  of  prolactin  [16]. 

Plasma  prolactin  Increase  Is  observed  2  min  after  Initiation  of  stress; 
15  to  20  minutes  after  the  stress,  the  concentration  of  plasma  prolactin  returns 
to  the  basal,  resting  level.  Prolactin  release  Is  also  under  the  influence  of 
catecholamine  levels  [12]. 

Although  the  functional  Importance  of  prolactin  release  remains  coscure 
(essential  actions  of  prolactin  are  mammotrophlc  and  lactogenic),  ’t  is  <ncwn 
that  this  nomcne  Is  released  during  stress  and  the  release  Is  mediated  by  the 
hypothalamus  [17]. 
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II.  MATERIALS  AND  tCTUOOS 


For  this  study#  the  concentration  of  plasma  prolactin  was  chosen  as  a 
sensitive  Indicator  of  possible  environmental  stresses  Induced  by  RFR.  To 
detect  and  quantitatively  evaluate  possible  Increases  In  plasma  prolactin  levels 
Induced  by  long-term  exposure  to  RFR»  blood  (0.3  mL)  was  periodically  sampled 
from  62  exposed  and  54  sham-exposed  animals.  Analysis  of  the  data  obtained  from 
blood  sample  assays  determined  whether  there  were  any  RFR-Induced  changes  In 
plasma  prolactin  concentration. 

Animal  s .  Male  Sprague-Oawl  ey  rats  were  used  In  this  study.  All 
experimental  animals  were  obtained  from  the  same  building  and  room  at  CAMM 
Research  Labs#  Wayne#  New  Jersey.  The  animals#  weighing  approximately  60  g, 
were  delivered  to  Emory  University  where  they  were  caged  singly  and  given  water 
and  food  (Purina  Rat  Chow)  ad,  libitum.  Temperature  In  the  animal  rooms  was 
maintained  at  24  +  1  °C  and  the  photo  period  was  12  hours/12  hours#  with  the 
lighted  phase  occurring  between  8  AM  and  8  PM. 

Experimental  Facility.  The  Georgia  Tech  Research  Institute's  Radiation 
Facility  [18]  consisted  of  8  collocated  rooms  on  the  basement  floor  of  the  Baker 
3u11d1ng  on  the  main  campus.  These  8  rooms  provided  a  closed#  complete  facility 
for  long-term  blceffects  studies  Involving  rodents. 

The  100  exposure  and  100  sham-exposure  animals  were  housed  In  two 
Identical#  collocated  rooms  In  the  Radiation  Facility.  Each  room  contained  a 
stack  of  circular#  parallel-plate  waveguides  fed  by  a  si otted-cyl inder  antenna 
system  for  radiating  the  animals.  The  stacks  of  parallel  waveguides  consisted 
of  five  3.6-m  (12  ft  )  diameter  plates  that  made  up  4  sets  of  circular 
waveguides.  Twenty-five  Individually  housed  rats  were  positioned  around  the 
c J rcumference  of  each  waveguide  set.  The  walls  of  both  rooms  were  I'ned  with 
anecholc  absorbing  material  and  shielded  with  aluminum  foil  to  prevent  excessive 
microwave  leakage  radiation. 

The  circular#  pa  rail  el -pi  ate  waveguide  assembly  provided  a  1.0  mW/cm- 
exposure  field  around  the  circumference  of  the  plates.  The  45.7-cm  (13  In.) 
slate  separation  distances  permitted  propagation  of  a  "E.q  mode  wave  Wth 
horizontal  polarization.  The  result  was  an  electric  field  vector  or'ented 
parallel  to  the  rat's  longitudinal  axis#  thereby  maximizing  the  coupling  between 
the  electric  field  and  the  rat.  The  power  density  displayed  a  ccs'ne-squared 
dependency  between  the  plates#  with  the  maximum  power  density  occurring  midway 
between  each  set  of  plates. 
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A  slotted-cyl inder  antenna  with  the  proper  diameter#  thickness#  slot 
length#  and  slot  width  dimensions  fed  the  stack  of  circular  waveguides  In  a 
manner  that  provided  an  essentially  constant  electric  field  Intensity  In  the 
azimuth  plane. 

Cages.  The  animal  cages  were  constructed  of  Plexiglas  to  facilitate  visual 
observation  of  the  rats  and  provide  radiofrequency  (RF)  transparency.  Each  cage 
was  22.9-cm  (9  In.)  long  by  12.7-cm  (5  In.)  wide  by  17.3-cm  (7  In.)  tall.  These 
dimensions  complied  with  recommended  caging  requirements  [18]  for  long-term 
housing  of  rats.  The  food  hopper  and  water  bottle  were  placed  on  the  distal 
side  of  the  cage  to  minimize  their  Interaction  with  the  exposure  field.  The 
glass  floor  rods  In  the  cage  were  oriented  perpendicular  to  the  cage's  long  axis 
to  encourage  the  rats  to  preferentially  align  themselves  parallel  to  the 
electric  field  vector,  Slpper  tubes  for  the  water  bottles  were  made  of  glass  to 
be  nonperturbing  In  the  field.  Evaluations  of  the  cages  conducted  In  the 
circular#  parallel-plate  waveguide  assembly  showed  field  scattering  from  the 
Plexiglas  and  water  to  be  below  the  range  of  detection. 

The  Radiation  Facility  also  contained  a  data  acquisition  system  for  storing 
and  processing  experimental  data#  an  electronic  balance  for  weighing  the  rats 
during  the  study#  and  rooms  for  transmitter  operation#  blood  sampling,  cage 
washing,  and  materials  storage. 

The  entire  Radiation  Facility  was  locked  to  avoid  unauthorized  entry.  This 
step  significantly  reduced  the  Introduction  of  noise  that  otherwise  could  have 
caused  artifacts  In  the  study  results.  Only  the  animal  caretaxer  and  the 
technician  who  sampled  blood  from  the  animals  were  permitted  uncontrolled  entry 
to  the  Facll Ity . 

Cannul atlon.  To  use  each  animal  as  Its  own  control#  arterial  blood  was 
sampled  by  means  of  Implanted  aortic  cannulas.  Cannulatlon  provided  a  simple# 
inexpensive  technique  that  permitted  remote#  stress-free  blood  sampling  In 
conscious#  unrestrained#  and  resting  rats  [2,203.  Arterial  blood  drawn  from  the 
chronically  Implanted  aortic  cannulas  was  assayed  for  plasma  prolactin.  Venous 
blood  was  not  sampled  because  the  blood  flow  in  veins  Is  laminar  and#  therefore# 
flews  In  discrete  layers  that  do  not  mix.  Only  physiologically  minute  amounts 
of  arterial  blood  (up  to  0.3  mL)  were  withdrawn  from  resting  rats  approximately 
every  2  weeks. 


The  carotid  artery  of  each  rat  was  cannulated  using  a  PE-10  cannula  8  to  10 
days  before  the  animals  entered  the  study.  The  surgery#  which  required  about  8 
min#  was  performed  using  ketamlne-xylazlne  anesthesia  (1:1  mixture;  ketamine  100 
mg/mL»  xylazlne  20  mg/mL>  l.m.  0.1  mL/100  g  of  body  weight).  The  cannulas  were 
filled  with  slightly  heparinized  saMne*  and  their  distal  ends  were  sealed  with 
nylon  plugs.  Stress  hormone  levels  returned  to  the  basal  values  about  3  days 
after  Implantation  of  the  chronic  arterial  cannulas  (Table  1).  The  first  blood 
sampling  occurred  10  days  after  aortic  cannulatlon. 

Blood  Sampling.  Restraint  and  handling  increase  stress  hormone  levels  In 
rats#  as  confirmed  during  the  study  (Table  2).  However#  the  animals  had  to  be 
handled  upon  removal  from  their  exposure  cage  and  placement  In  the  "sampling 
box"  In  preparation  for  blood  withdrawal.  To  avoid  the  undesired  effects  of 
handling  and  stress  on  hormone  levels#  blood  from  the  aortic  cannula  was  sampled 
30  min  after  the  animal  was  placed  In  the  sampling  box.  This  procedure 
permitted  the  altered  plasma  prolactin  level  sufficient  time  to  return  to  Its 
basal  value  (Table  2).  Each  animal  was  previously  preconditioned  for  the 
sampling  box  through  a  regime  of  several  30-mln-long  experiments  conducted 
during  a  1-week  period  before  entering  the  study. 

After  acclimating  for  30  min  In  the  sampling  box#  the  rat's  cannula  was 
positioned  through  the  slot  In  the  top  of  the  box  (Fig.  1).  The  heparinized 
saline  was  then  removed  from  the  cannula#  and  a  0.3-mL  blood  sample  was  taken 
from  the  resting  rat.  The  withdrawal  of  larger  amounts  of  blood  from  the 
cannulated  rats  would  have  altered  the  level  of  stress  hormones.  Using  a 
sterile  1-cm3 tubercul In  syringe  fitted  with  a  30-ga  needle#  the  blood  sample  was 
taken  from  the  cannula.  The  syringe  and  the  needle  were  rinsed  with 
ethyl enedl ami netetraacetate  (EDTA)  before  sampling.  The  blood  sample  was  placed 
In  an  EDTA-treated  0.3-mL  capillary  blood  collection  container  (Walter  Sarstedt 
Co.»  Princeton,  New  Jersey)#  shaken#  and  then  placed  on  ice.  The  blood  sampling 
procedure  required  about  2  min  for  each  rat. 

Blood  Sampling  Schedule.  Figure  2  shows  the  sampling  schedule  designed  for 
the  experiment.  Note  that  the  200  rats  were  Introduced  into  the  study  in  4 
groups  of  50  animals  each.  The  groups  entered  In  a  staggered  manner  to 
facilitate  the  process  of  logging  In  and  establishing  the  new  animals.  Each 
group  contained  25  exposure  and  25  sham-exposure  animals.  Of  the  25  exposure 


*  0.5  cm3  heparin  sodium  (from  beef  lung)#  1000  units/mL  per  30  cm3  sal  1 ne. 
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TABLE  1 


TABLE  2 


PLASMA  PROLACTIN  VALUES  (ng/mL)  ±$0 
OAYS  AFTER  IMPLANTATION  OF  THE 
SAMPLING  (0.3  mL)* 

Arterial  Blood  Sampled  After  Days 

OBTAINED  IN  RESTING  RATS  SEVERAL 
CHRONIC  AORTIC  CANNULA  FOR  BLOOD 

of  Aortic  Cannulatlon 

I.  dav 

3  davs 

5  davs 

7  davs 

14  davs 

58 

16 

4 

11 

30 

32 

21 

23 

22 

7 

47 

21 

15 

12 

7 

30 

7 

12 

6 

12 

12 

19 

19 

21 

10 

52 

46 

14 

13 

18 

38 

29 

16 

10 

9 

18 

41 

9 

3 

15 

27 

7 

7 

18 

17 

42 

_a 

21 

_2 

-2 

36  +  15 

22  +  14 

14  +  6 

12  +  7 

13  +  7 

•Each  group 

of  rats  consisted  of  10  animals 

(1>3>5>7»  and 

14  days). 

The  animals 

were  adapted  to 

their  cages 

for  3  weeks 

before  the 

cannulas  were  Implanted. 

PLASMA 

PROLACTIN  VALUES  (ng/mL) 

+SD 

IN  10 

RESTING  RATS  AND  IN  THE  SAME 

RATS  7 

MIN  AFTER  PLACEMENT  INTO 

NEW 

CAGES* 

7 

min  After  Placement 

1 

25 

29 

2 

12 

21 

3 

8 

16 

4 

11 

43 

5 

14 

29 

6 

12 

11 

7 

17 

38 

8 

18 

29 

9 

7 

26 

10 

12 

21 

14  +  5 

27  +  10 

•Sampling  through  chronic  aortic  cannula  while  the  animal  rests  In  Its 
home  cage. 


Figure  1.  Sampling  0.3  mL  of  blood  from  the  chronically 
implanted  aortic  cannula  of  a  resting  rat. 


(or  sham-exposure)  animals,  20  were  sampled  for  plasma  stress  hormones,  while 
the  remaining  5  were  used  for  hematology  studies. 

The  sampling  duration  was  36  weeks  long,  including  a  6-week  preexposure 
adaptation  period,  a  24-week  exposure  period,  and  a  6-week  postexposure  period. 
With  allowing  for  group  staggering,  the  experiment  duration  was  42  weeks  long 
(since  the  4  groups  entered  2  weeks  apart  from  one  another).  Plasma  prolactin 
was  sampled  for  all  periods  marked  (A)  in  Figure  2.  Therefore,  each  animal 
should  have  been  sampled  for  plasma  prolactin  at  weeks  -6,  -3,  0,  3,  6,  ...,  27. 
This  schedule  was  rather  rigorous,  and  therefore  could  tolerate  slight 
fluctuations  in  protocol  without  ill  effects. 

Prolactin  Determination.  Rapid,  sensitive,  and  specific  radioimmunoassays 
that  require  a  minimum  quantity  of  blood  were  used  In  this  study.  These 
qualities  were  especially  important  because  repetitive  sampling  was  required  and 
small  laboratory  animals  were  used.  Plasma  prolactin  from  individual  plasma 
samples  was  measured  In  duplicate  by  double  antibody  radioimmunoassays  for  rat 
prolactin  [21]  using  the  NIAMD*  reagents.  Results  are  expressed  as  ng/mL,  and 
the  reference  standard  was  Rat  Prolactin  RP-2. 

Resting  Value  of  Plasma  Prolactin.  At  the  initiation  of  the  study, 
preliminary  experiments  were  performed  to  determine  the  basal  value  of  plasma 
prolactin  In  the  cannulated  resting  rats  (see  Table  1).  Plasma  stress  hormones 
In  the  rat  follow  circadian  rhythm.  Increasing  during  evening  hours  and 
decreasing  to  the  lowest  level  between  9  AM  and  1  PM  [22,23].  To  avoid  the 
effects  of  circadian  rhythm  on  the  study  results,  blood  sampling  occurred  only 
between  9  AM  and  1  PM  when  prolactin  concentration  was  at  its  lowest  (true 
resting)  level  [9]. 
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III.  RESULTS  AND  ANALYSIS 

Appendix  A  contains  the  data  collected  during  the  course  of  the  pre¬ 
exposure  and  radiation  periods  for  both  exposure  and  sham-exposure  animals. 
Over  the  entire  blood  sampling  period#  there  was  considerable  variance  In  the 
data#  suggesting  animal  activity  at  the  time  of  blood  sampling.  Since  the 
sampling  boxes  had  opaque  walls#  the  physical  activity  of  each  animal 
Immediately  prior  to  sampling  was  not  recorded;  however#  each  animal  had 
sufficient  time  (30  min  or  more)  to  return  to  basal  hormonal  level  after  the 
stimulation  of  being  placed  Into  the  sampling  box. 

Figures  3  and  4  present  the  data  of  Appendix  A  In  scatter  form  (one  plot 
corresponds  to  sham-exposure  animals#  the  other  to  exposure  animals).  The 
dotted  line  passes  through  the  mean  hormone  response  at  each  week.  In  general, 
plots  of  both  exposure  and  sham-exposure  hormone  concentration  versus  time  were 
essentially  linear  (although  there  was  some  curvature  present  at  the  exposure 
onset  and  conclusion).  Furthermore,  the  trend  of  the  data  suggested  that  plasma 
prolactin  concentrations  In  both  exposure  and  sham-exposure  groups  began 
somewhat  high,  declined  Into  the  study#  and  then  rose  slightly  toward  the  end  of 
the  exposure.  There  was  little  variation  In  the  two  plots  when  they  were 
overlaid  and  compared  (Fig.  5).  This  was  preliminary  evidence  Indicating  that 
435-MHz  RFR  did  not  Increase  resting  plasma  prolactin  concentrations.  To  attach 
numerical  probabilities  to  this  conclusion#  the  data  were  statistically 
anal yzed. 

The  plasma  prolactin  data  were  analyzed  with  linear  regression  model- 
building  techniques.  A  quadratic  model  (hormone  concentration  as  a  function  of 
time)  was  constructed  to  fit  the  data.  Terms  of  the  quadratic  model  were  then 
tested  to  determine  whether  or  not  there  were  significant  microwave-induced 
effects  on  hormone  concentration.  Appendix  B  contains  a  detailed  discussion  of 
the  methodology#  procedure#  and  results  of  the  statistical  analysis. 

Results  of  the  analysis  Indicated  that#  If  there  were  any  RFR-Induced 
effects  on  plasma  prolactin  concentration#  these  effects  were  within  the  range 
of  +  3.32  ng/mL  from  the  estimated  normal  resting  value  of  17.05  ng/mL.  Since 
this  range  was  within  the  normal  range  of  plasma  prolactin  concentration 
variability  In  unstressed  male  rats#  there  was#  from  a  practical  standpoint#  no 
Indication  of  RFR-Induced  stress  affecting  animal  resting  plasma  prolactin 
concentrations. 
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Time  (in  weeks) 

NOTE:  3677  OBS  HAD  MISSING  VALUES  OR  WERE  OUT  Of  RANGE  2A5  OBS  HIDOEN 

Figure  3.  Prolactin  concentration  data  scatter  diagram  (sham-exposure  group). 
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Time  ( in  weeks) 

3780  OBS  MAD  MISSING  VALUES  OR  WERE  OUT  OF  RANGE  248  OBS  HIDDEN 

Figure  4.  Prolactin  concentration  data  scatter  diagram  (exposure  group). 
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Figure  5.  Mean  plasma  prolactin  concentrations  versus  time. 


IV.  DISCUSSION 


It  Is  known  that  stress  increases  the  level  of  plasma  stress  hormones. 
Thus#  handling  of  the  animals  C24]#  exercise  [25]#  Immobilization  [26], 
withdrawal  of  large  volumes  of  blood  [27]#  exposure  to  new  or  unfamiliar  housing 
[28]#  noise#  hypoxia  [29]#  cold  or  heat  exposure  [30]#  and  many  other 
environmental  factors  Increase  the  plasma  concentration  of  stress  hormones. 
Both  neurogenic  (emotional)  stimuli  and  systemic  (somatic)  stimuli  are  effective 
In  evoking  Increased  secretion  of  stress  hormones  in  animals  (and  In  man)#  and 
these  stimuli  had  to  be  avoided  In  the  Radiation  Facility  used  for  this  study. 

Handling  and  removing  a  rat  from  Its  cage  also  Induces  an  Increase  In 
plasma  stress  hormones  even  If  the  stimulus  is  removed  Immediately.  The 
Increase  In  plasma  stress  hormones  was  observed  for  20  to  30  min  [31]. 

The  plasma  prolactin  Increase  In  response  to  stress  [4,5,12]  can  be 
quantified.  The  degree  of  plasma  prolactin  Increase  Is  related  to  the  type  and 
Intensity  of  stress  to  which  the  animal  Is  exposed  [32]  as  well  as  to  the 
duration  of  stress.  The  physiological  Importance  of  increased  prolactin  release 
In  response  to  stress  remains  poorly  understood. 

There  are  few  studies  that  deal  with  the  effects  of  long-lasting  stress. 
Burchfield  et  al  .  [33]  demonstrated  that  the  resting  plasma  corticosterone 

levels  in  chronically  stressed  rats  had  elevated  as  much  as  in  control  animals 
during  acute  stress#  but  plasma  ACTH  levels  remained  unchanged.  In  another 
study#  it  was  shown  that  adaptation  to  stress  did  not  result  In  an  increased 
rate  of  adrenocortical  response  and  "an  overall  Increased  responsiveness  cf  the 
pituitary-adrenal  system”  [34]. 

The  high  sensitivity  of  the  braln-pituitary-prolactln  system  observed 
during  stress  demands  that  blood  sampling  be  done  remotely.  Repeated  sampling 
of  blood  from  the  same  cannulated  rat  provided  reliable  resting  patterns  of 
prolactin  secretion  that  would  reveal  any  increases  induced  by  a  long-term#  Icw- 
level  RFR  environment.  Apparently#  even  the  smallest  environmental 
perturbation#  such  as  low-level  RFR,  would  be  detec  table  if  It  had  any 
significant  Influence  on  the  release  of  this  hormone.  Although  relocation  of  a 
rat  from  the  cage  Into  the  sampling  box  30  min  before  blood  sampling  slightly 
disturbed  the  environment  of  the  rat#  such  perturbations  did  not  alter  resting 
plasma  prolactin  levels  at  the  time  of  sampling  (Table  2). 
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Results  of  a  study  concerning  plasma  ACTH  and  corticosterone  concentrations 
In  rats  exposed  for  6  months  to  the  same  RFR  environment  used  In  this  study  were 
reported  previously  [2].  These  results  showed  that  plasma  corticosterone  and 
plasma  ACTH  concentrations  were  not  changed  In  rats  exposed  to  low-level,  pulsed 
RFR  fields  for  a  6-month  duration.  In  this  report,  plasma  prolactin  levels  In 
the  same  animals  are  reported.  These  3  hormones  were  studied  because  every 
stress  does  not  release  all  stress  hormones.  While  In  certain  cases, 
associations  are  observed  In  the  release  of  some  stress  hormones  (for  Instance, 
corticosterone  and  ACTH,  [35]),  multiple  hormone  release  Is  not  always  observed. 
Furthermore,  while  corticosterone  Is  released  In  a  pulsatile  fashion,  the 
release  of  prolactin,  though  also  pulsatile,  induces  smaller  variations  from  the 
mean  and  thus  might  provide  a  better  method  for  measuring  the  resting  level  of 
stress  hormones. 

It  has  already  been  shown  that  short-term  exposure  to  low-level  microwave 
radiation  does  not  change  the  plasma  level  of  some  stress  hormones  In  rats 
[36,37].  Johnson  and  associates  [38]  found  an  elevation  of  plasma 
corticosterone  the  first  time  the  blood  was  sampled  from  microwave-exposed  rats 
In  their  long-term  study.  In  the  same  study,  plasma  corticosterone  returned  to 
resting  control  levels  throughout  the  remaining  2-year  period. 

As  previously  mentioned,  plasma  prolactin  Is  a  sensitive  Indicator  of 
various  types  of  environmental  stress  in  mammalian  systems.  Stressors  lead  to 
Increased  prolactin  release  and  an  Increased  plasma  prolactin  concentration. 
This  Increase  depends  on  the  Intensity  and  duration  of  stress,  and  can  reach  8 
to  10  times  the  normal  resting  plasma  concentration.  Our  results  show  that  Icw- 
level  RFR  does  not  change  plasma  prolactin  levels  In  rats.  The  statistical 
analysis  Indicates  that  any  RFR-Induced  effects  on  rat  resting  prolactin 
concentration  would  lay  within  a  range  of  +3.32  ng/ml  from  an  estimated  resting 
concentration  of  17.05  ng/mL.  These  values  are  not  typical  of  rats  exposed  to 
stress.  Therefore,  this  study  concludes  that  a  1.0  mW/cm^  435-MHz  pulsed-wave 
(1.0  us  width,  1  kHz  pulse  rate)  RFR  environment  did  not  Induce  any  detectable 
Increase  In  stress,  as  measured  by  resting  prolactin  concentration,  in  the 
exposure  group  of  62  cannulated  male  Sprague-Oawl ey  rats  when  compared  to  a 
sham-exposure  group  of  64  cannulated  male  Sprague-Cawley  rats. 


Y.  REFERENCES 


Toler,  J.»  Bonasera,  S.»  and  Popov1c,V.  Long-term  Bioeffects  of  435-MHz 
Radiofrequency  Radiation  on  Selected  Blood-borne  Endpoints  In  Cannulated 
Rats.  Volume  1.  Engineering  considerations.  USAFSAM-TR-86-4,  September 
1986. 

PopovIc.V.,  Toler, J.,  Bonasera, S.,  PopovIc.P.,  Honeycutt, C. #  and  Sgoutas.D. 
Long-term  Bioeffects  of  435-MHz  Radiofrequency  Radiation  on  Selected  Blood- 
borne  Endpoints  In  Cannulated  Rats.  Volume  2.  Plasma  ACTH  and  plasma 
corticosterone  results.  USAFSAM-TR-37-5  (in  preparation  ). 

Freeman,  M.E.,  and  Nell.J.D.  The  Pattern  of  Prolactin  Secretion  During 
Pseudopregnancy  In  the  Rats  A  Dally  Nocturnal  Surge.  Endocrinology 
90:1292-1294  (1972). 

Matthel J,  J .A. ,  and  Pljkeren,  T.A.  Plasma  Prolactin  In  Undisturbed 
Cannulated  Male  Rats:  Effects  of  Perphenazine  Frequent  Sampling  Stress  and 
Castration  Plus  Oestrone  Treatment.  Acta  Endocrlnologlca  84:51-61  (1977). 

Plercy.M.,  and  Shin,  S.H.  Comparative  Studies  of  Prolactin  Secretion  In 
Estradiol-Primed  and  Normal  Male  Rats  Induced  by  Ether  Stress,  Pimozide  and 
TRH.  Neuroendocrinology  31:270-275  (1980). 

Stern,  J.M.,  Thomas,  D.  A.,  Rabll.J.,  and  Barf lei d ,  R.J  .  Do  Pup  Ultrasound 
Cries  Provoke  Prolactin  Secretion  In  Lactatlng  Rats?  Hormones  and  Behavior 
18:86-94  (1984). 

Fromm,  M.,  Oelkers,  W.»  and  Hegel,  U.  Time  course  of  aldosterone  and 
corticosterone  plasma  levels  in  rats  during  general  anaesthesia  and 
abdominal  surgery.  Pflugers  Arch  399:249-254  (1983). 

Hlroshlge.T.  Hormonal  rhythm  and  feeding  behavior.  J  Autonom  Nerv  System 
10:337-346  (1984). 

Matthel j, J .A. M. ,  and  Swarts, H. J .M.  Circadian  Variations  of  the  Response  of 
Plasma  Prolactin  to  Brief  Exposure  to  Ether  In  the  Male  Rat.  Acta 
Endocrlnologlca  94:327-331  (1980). 

Deeter,  W.T.,  and  Mueller,  G.P.  Differential  Effects  of  Warm-  and  Cold- 
Ambient  Temperature  on  81ood  Levels  of  2 -Endorphin  and  Prolactin  In  the 
Rat.  Soc  Exper  Biol  Med  168:369-372  (1981). 

Fayez.M.,  Ahmed,  H.H.,  and  Quadrl,  S.K.  Effect  of  Proplcmazlne  on  Plasma 
Prolactin  in  the  Rat,  Counteraction  by  L-Dopa.  Proc  Soc  Exper  Biol  Med 
178:606-609  (1985). 

Tobias, H.,  Carr.L.,  and  Voogt,J.  Catecholamine  Mechanisms  In  the  Feedback 
Effects  of  Estradiol  Benzoate  on  the  Release  of  LH  and  Prolactin.  Proc  Soc 
Exper  Biol  Med  174:284-290  (1983). 

OkaJIma.T.,  Motomatsu,  T. ,  Kato.K.,  and  Ibayashl,  H.  The  Stimulatory  Effect 
of  Beta-Endorphin  on  the  Plasma  Prolactin  Levels  was  Diminished  In  the  Rats 
Treated  with  6-Hvdroxydooamlne.  Life  Sciences  26:699-705  fl98m. 


trwsTrrwsrjrvTrurir 


w.iwwfw^injtMwwwwwwniwumirjwjwii'w.'wwiwwpjwiiyjyjrjif.'Hjpjr.'T.'P.'ww.'PJWW^ 


I 

J 


i 

I 

14.  Siegel,  R. A.*  Confortl.N.,  and  Chowers, I.  Neural  Pathways  Mediating  the  *J 

Prolactin  Secretory  Response  to  Acute  Neurogenic  Stress  In  the  Male  Rat.  .1 

Brain  Research  198:43-53  (1980).  N 

15.  Col ombo, J .A. ,  and  Luth,  P.E.  Immediate,  Preferential  Prolactin  Release  j] 

After  Discrete  Brain  Lesions  In  Male  Rats.  Exper  Neurol  72,  637-697  I 

(1981).  3 

16.  Seggle,  J.  Amygdala  Lesions  and  24-h  Variation  In  Plasma  Corticosterone, 

Growth  Hormone,  and  Prolactin  Levels.  Can  J  Physiol  Pharmacol  53:249-253  'l 

(1980).  i 


17.  Mloduszewskl,  R.,  and  Critchlow,  V.  Marked  and  Selective  Suppression  of 

Prolactin  Secretion  After  Interruption  of  Neural  Connections  of  the  Medial 
Basal  Hypothalamus  of  the  Rat.  Brain  Research  216:35-44  (1981). 

18.  OHE'W  Publication  (NIH)  78-  23.  Guide  for  the  care  and  use  of  laboratory 
animals.  Revised  1978. 

19.  Toler.J.,  and  Popovic.V.  Operational  evaluation  of  a  new  435— MHz  raciation 
facility.  USAFSAM-7R-33-21,  October  1983. 

20.  Popovlc,  V.,  and  Kent,  K.M.  120-day  study  of  cardiac  output  In 

unanesthetized  rats.  Am  J  Physiol  207:767-770  (1964). 

21.  Nlswender,  G.D.,  Chen,  C.L.,  Mldgely,  A.R.,  Meltes,  J.,  and  Ellis,  J. 
Radioimmunoassay  for  Rat  Prolactin,  Proc  Soc  Exp  Biol  yea  130:793-797 
(1969). 

22.  Guillemin,  R.,  Dean,  W.E.,  and  Lieoelt,  R.A.  Nychthemeral  variations  in 

dlasma-fr9e  corticostercid  levels  of  the  rat.  Prcc  See  Exp  Biol  yeo 
101:394-395  (1959). 


! 

/" 

/ 

r 

/ 

/ 

/ 

\ 


23.  Simcn,  M.L.,  and  George, R.  Ciurnal  variations  In  plasma  corticosterone  anc 

growth  hormone  as  correlated  with  regional  variations  in  norep  1 neDhri ne, 
dooamine  and  serotonin  contents  of  rat  brain.  Neuroendocrinology  17:125- 
133  (1975). 

24.  <n1gge,  K.M.,  Penrod,  C.H.,  and  Scninoler,J.  In  vitro  and  v'vo  adrenal 

corticosterone  secretion  following  stress.  Am  J  Physiol  196:579-532 

(1959). 

25.  Hennessy,  J .W. ,  Levin,  R.,  and  Levine, S.  Influence  of  riential  factors  anc 
gonadal  hormones  on  pituitary-adrenal  responses  of  the  mouse  to  novelty  and 
electric  shock .  J  ^ysiol  Dsychol  51:773-7'”  '197''. 

26.  Baron, S.,  and  Brush, F.R.  Effects  of  acute  and  chronic  restraint  and  estrus 
cycle  on  pituitary-adrenal  function  In  the  rat.  Hormones  Senav  12:213-224 
(1979). 

27.  Barton,  R.N.,  and  Passingham,  3.J .  Early  resconses  to  hemorrhage  In  the 

conscious  rat:  Effects  of  corticosterone.  Am  J  Physiol  (Regulatory 

Integrative  Comp  Physiol  12)  243 :R415-R423,  (1982). 


■> 

h 

,*< 

r> 

h 

* 

hj 

* 


•J 


.h 


28.  Fagln, K.D. *  and  Neill,  J.D.  Involvement  of  the  neuro Intermediate  lobe  of 

the  pituitary  gland  In  the  secretion  of  prolactin  and  luteinizing  hormone 
In  the  rat.  Life  Scl  30:1135-1141  (1982). 

29.  0u,  L.C.»  and  Tenney,  S.M.  Adrenocortical  function  In  rats  chronically 

exposed  to  high  altitude.  J  Appl  Physiol:  Resplrat  Environ  Exercise 

Physiol  47:1185-1187  (1979). 

30.  Gale.G.C.#  Hayward, J .S. #  Green, W.L.,  Wu.S.Y.,  Schll ler, H.J . ,  and  Jackson,!. 

Endocrine  responses  to  acute  cold  In  man.  Fed  Proc  34:301  (1975). 

31.  Seggle,  J  .A. ,  and  Brown,  G.M.  Stress  responses  patterns  of  plasma 

corticosterone#  prolactin#  and  growth  hormone  In  the  rat#  following 
handling  or  exposure  to  novel  environment.  Can  J  Physiol  Pharmacol  53:629- 
637  (1975). 

32.  Bassett, J.R.#  Cal rncross,K .0. ,  and  King, M.G.  Parameters  of  novelty,  shock 

predictability  and  response  contingency  to  corticosterone  release  In  the 
rat.  Physiol  Behav  10:901-907  (1973). 

33.  Burchfield,  S.R.,  Woods,  S.C.,  and  Ellch.M.S.  Pituitary  adrenocortical 
response  to  chronic  Intermittent  stress.  Physiol  Behav  24(2) :297-302 
(1980). 

34.  Sakell  arls,  P.C. #  and  Vernlkos-Oannell Is,  J.  Increased  response  of  the 

pituitary-adrenal  system  in  rats  adapted  to  chronic  stress.  Endocrinology 
97:597-602  (1975). 

35.  0oe11,  R.G.#  Dallman.M.F.,  Cl ayton ,  R.B. ,  Gray,  G.D.,  and  Levine,  S. 

Dissociation  of  adrenal  corticosteroid  production  from  AC^H  In  water- 
restricted  female  rats.  Am  J  Physiol  241  (Regulatory  Integrative  Comp 
Physiol  10):  R21-R24  (1981). 

36.  Lovely,  R.H.,  Myers,  D.E.,  and  Guy,  A.W.  Irradiation  of  rats  oy  918  '•'Hz 

microwaves  of  2.5  mW/crrr:  Delineating  the  dose-response  rel aticnsni p. 

Radio  Scl  12:139-146  (1977). 

37.  Johnson,  R.3.#  Mlzumorl,  S.#  Lovely,  R.H.#  and  Guy,  A.W.  Adaptation  to 
microwave  exposure  as  a  function  of  power  density  and  ambient  temperatures 
In  the  rat.  In  electromagnetic  fields  In  biological  systems.  Abstracts 
1978  IMPI/URSI  Symposium,  Ottawa#  Canada, 27-30  June  1978. 

38.  Johnson, R.B. #  Spackman.D.#  Crowley, J.,  Thompson, D.#  Chou.C.K.#  Kunz.L.L.# 
and  Guy,  A.W.  Effects  of  long-term  low-level  radiofrequency  radiation 
exposure  on  rats.  Vol .  4.  Open-field  behavior  and  corticosterone.  5eoort 
USAFSAM-TR-83-42 ,  December  1983. 


17 


i 


Prolactin 
Control  II 


!—■■■  ■■■■■■■■■■■ ■■mal 

Bmwbbbbwbmmbmhb 

— — BIMH—M— — I 


Sf)  17' U  '2  I 


■I  7, 


to  ^  to  io 


J4  70'I3Z 


70 


Wto'b'to'ty 


Tffi 


qu'W  xik 


nm 


i  m 

me 


SBsmm 


9 


SSm 


r?ai—r—EiHiuF»-j 


!— — E— — KMIHW— J 


I— ■ETiBME— — — BH— — — 

!■— — — — — —Sami 

I— — — — MHI— ■■■ay 

— ■ —Hp — — he— — —a— —I 

I—  — a— — ——mm 

■—  — — — — ■— ■— —aa ■ 

II—  E— — — MIBjllJ 


Prolactin 
MW  IV 


IM  t  trmm  ■**  -KK  KKUKKKUKMKttKKKKKKKKK  UK  UIK  1»1M  1UK  IKK  IKK  IKK  IKK  i«K  »K  IHK  IKK  1  »K  W  .1  .1 


APPOOIX  B 

STATISTICAL  >CTHOOOLOGY 


! 

t 

i 


25 


'  *  «v"  ^  A 


APPENDIX  B 


STATISTICAL  JCTHOOOLOGY 

The  balanced  design  of  this  experiment  (requiring  that  25  animals  from  each 
group  be  sampled  once  every  3  weeks  for  stress  hormones)  should  have  produced 
data  easily  tested  by  balanced*  2-way  analysis  of  variance  (ANOVA)  statistics 
with  12  levels  of  factor  A  (time)  and  2  levels  of  factor  B  (RF  radiation). 
However*  data  collection  did  not  proceed  according  to  protocol  In  that*  In 
numerous  cases*  samples  were  collected  at  odd  Intervals  (Invalidating  the 
orthogonality  of  the  design)  and  the  number  of  samples  taken  per  week  varied 
more  or  less  than  25  (unbalancing  the  design).  These  two  factors  combined  to 
lower  the  power  of  ANOVA  statistics  (power  Is  defined  as  the  ability  to  reject 
the  null  hypothesis  given  the  null  hypothesis  should  be  rejected)  trying  to  test 
the  model 

»Uk  -  "  ♦  T1  *3j  +  rSu  *  £<Jk'  ,B‘1) 

hormone  concentration  (response), 
the  normal  hormone  resting  concentration, 
the  change  In  hormone  resting  concentration  Induced  by  RFR, 

the  change  In  hormone  resting  concentration  Induced  by  time, 

the  change  In  hormone  resting  concentration  induced  by  the 

Interaction  between  RFR  and  time,  and 
noise  within  the  system  (sampling  and  assaying  errors) 

for  the  following  hypotheses: 

H0:  T  0  =  T  1  =  °» 

t g  or  0  (RFR-Induced  effects),  (3-2) 

H0:  3  i  =  32  =  •  •  •  ■  3  12  =  °# 

H^:  at  least  one  3  j  +0  (time-induced  effects),  (B— 3 ) 

Hg:  t3jj  =  0,  and 

H^:  at  least  one  t3^^  0  (interaction  between  RFR  and  time).  (B-4) 
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However*  examination  of  the  collected  data  suggested  an  alternative 
approach  In  that  the  data  resembled  what  might  have  been  collected  In  an 
unplanned  experiment  monitoring  over  time  the  operation  (In  this  case, 
characterized  by  resting  animal  hormone  concentrations)  of  an  established  RF 
radiation  facility.  Oata  of  this  type  are  often  successfully  treated  by 
employing  linear  regression  techniques  to  develop,  build,  and  test  a  linear  (or 
Intrinsically  linear)  model  whose  parameters  can  be  used  to  predict  the  system 
response  at  various  treatment  levels.  Therefore,  we  decided  to  proceed  with  a 
regression  approach  to  data  analysis. 

The  first  step  In  the  regression  approach  to  data  analysis  was  to  define  an 
Initial  model  to  fit  the  data,  and  to  test  the  properties  of  this  model. 

Visual  Inspection  of  the  scatter  diagrams  of  Figures  3  and  4  showed  an 
essentially  linear  plasma  prolactin  response  versus  time.  Therefore,  there  was 
a  nonzero  3g  In  the  final  model,  and  tests  were  conducted  for  a  RFR-Induced 
effect  on  this  Intercept  with  the  term  aQz*  Also,  there  was  sufficient 
curvature  In  the  plot  (particularly  at  exposure  onset  and  termination)  to 
justify  the  Inclusion  of  linear  terms  (S^  and  a  ^z)  and  quadratic  terms  (3^  and 

%LlZ). 

The  initial  model  therefore  became: 

y  =  3o  +  5ix  +  3  n*2  +  3  oz  +  3izx  +  3nzx2  {B_5) 

where  y  -  the  plasma  prolactin  concentration, 
x  *  the  time  (in  weeks),  and 

z  *  a  categorical  variable  with  value  of  0  for  animals  In  the  sham- 
exposure  group  and  1  for  animals  In  exposure  group. 

At  this  point,  raw  data  from  the  prolactin  spreadsheet  were  put  on  computer 
file  (see  Appendix  A).  A  Statistical  Analysis  System  (SAS)  program  (see 
Appendix  C)  was  then  written  to  read  the  raw  data  file,  format  the  data  for 
analysis,  and  perform  a  variety  of  statistical  tests  on  the  model. 

The  first  test  Identified  terms  In  the  general  model  which  contributed  the 
least  to  forming  a  statistically  significant  regression.  Two  stepwise 
regression  procedures  were  employed:  forward  regression  and  maximum  R^ 

regression.  Forward  regression  procedures  entered  variables  Into  the  model  in 
such  a  way  as  to  produce  the  greatest  Increase  In  R^  (r2  t,eing  a  measure  of  the 


percentage  of  variation  In  the  data  set  which  Is  explained  by  the  statistical 
model)  while  ensuring  that  the  variable  entered  was  statistically  significant  at 
a  significance  of  0.15.  The  forward  stepwise  regression  produced  the  model  (see 
Appendix  0  for  the  SAS  forward  and  maximum  R2  analysis  output): 


+  S  .X  +3  ,  ,x  » 


(8-6) 


where  all  variables  were  as  previously  defined. 

2 

The  second  stepwise  procedure  employed  was  maximum  R*  regression  (MAXR). 
Maximum  R*  regression  functioned  essentially  the  same  way  as  the  forward 
procedure,  the  distinction  being  that  MAXR  entered  a  variable  Into  the  model  so 
long  as  the  Introduction  of  that  variable  Increased  the  R2  ratio  (even  If  the 
variable  was  found  to  be  otherwise  statistically  Insignificant).  Thus,  MAXR 
first  found  the  best  possible  2-parameter  model,  then  the  best  possible  3- 
parameter  model,  up  to  the  best  all -factor  model. 

The  combined  output  of  these  2  programs  gave  a  good  indication  (when  viewed 
with  estimates  of  the  coefficients  In  the  all-parameter  model)  of  which  terms  In 
the  original  model  could  be  removed  without  compromising  the  final  model’s 
predictive  power.  Both  forward  and  maximum  R2  regression  determined  that,  at 
the  0.15  significance  level,  neither  3q,  i  or  3^  were  important  to  the 
original  model.  (ig  was  significant  at  J  8  0.3818,  3^  was  significant  at  2 = 
0.8850,  and  3^  was  significant  at  i=  0.7264.)  Thus,  since  the  terms  modeling 
the  RFR  interaction  effect  were  Insignificant,  the  conclusion  was  drawn  that  RFR 
exposure  did  not  produce  a  detectable  effect  on  plasma  prolactin  concentrations. 

Note  that  the  estimated  values  for  3  g,  3  and  3^  were  all  found  to  be 
significant  at  a  level  greater  than  3  =  0.001  (Appendix  E,  page  44).  This 
Indicated  that  the  plasma  prolactin  concentrations  In  both  exposure  and  sham- 
exposure  groups  varied  over  the  duration  of  the  experiment.  This  curvature  from 
the  straight  line  case  (y  =  3g,  which  would  indicate  that  all  hormone 
concentrations  remained  constant  ever  time)  took  into  account  the  slightly 
higher  values  of  plasma  prolactin  at  the  experiment  onset  and  conclusion  (17  to 
19  ng/ml)  as  compared  to  plasma  prolactin  concentrations  in  the  middle  of  the 
study  (approximately  11  ng/ml).  These  predictions  are  rough  estimations  from 
the  model  since  the  confidence  Intervals  (provided  under  a  separate  cover)  on 
plasma  prolactin  concentration  were  the  same  width  as  normal  hormone  ranges  (10 
to  15  ng/ml). 


To  complete  the  analysis  (with  regards  to  the  question  of  RFR-Induced 
bioeffects)  required  the  subsequent  determination  of  the  maximum  perturbation  In 
resting  prolactin  levels  that  the  experimental  protocol  was  capable  of 
detecting.  However,  In  order  for  results  from  the  linear  regression  to  be 
considered  significant,  it  was  first  Important  to  verify  that  the  assumptions 
made  In  forming  the  linear  model  were  not  violated  during  the  model-building 
procedure.  These  assumptions  Included  no  lack-of-flt  In  the  model,  and  that  the 
residuals  from  the  fitted  model  followed  a  normal.  Independent  distribution 
(termed  NID  (0,  a^)).  First,  a  lack-of-flt  test  was  performed  on  the  data  by 
obtaining  (In  the  revised  model)  sum-of-squares  regression  error  and  sum-of- 
squares  pure  error.  Since  there  were  repeated  measurements  taken  at  each  week 
for  both  the  exposed  and  sham-exposed  animals.  It  was  therefore  possible  to 
break  the  model  sum-of-squares  error  Into  lack-of-flt  and  pure  error  terms. 
First,  the  model  sum-of-squares  error  was  obtained  by  running  a  regression  on 
the  revised  model  and  reading  the  term  from  the  resulting  ANOVA  table.  To 
obtain  a  sum-of-squares  pure  error  term,  the  SAS  General  Linear  Models  (GLM) 
procedure  was  applied  to  the  data  (33  levels  of  time  treatment,  2  levels  of  RFR 
treatment).  The  sum-of-squares  error  term  yielded  by  the  GLM  represented  a  sum- 
of-squares  pure  error  (due  to  sampling  variation)  In  the  regression.  Sum-of- 
squares  lack-of-flt  was  then  the  regression  sum-of-squares  error  minus  the  sum- 


of-squares  pure  error.  Calculations  to  compute  the  critical  value  Fg  from  these 


sum-of-squares  terms  are  detailed  in  Appendix  E. 


The  computed  test  statistic  Fg  exceeded  the  critical  value,  thereby 


indicating  significant  lack-of-fit.  Normally,  this  result  would  be  faintly 

disturbing  since  It  would  require  refitting  the  model  using  transformed  rather 

than  raw  data  values.  The  transformation  of  the  dependent  variable  y  was 

definitely  undesirable,  since  the  residual  plots  indicated  that  the  residuals  of 

2 

y  (using  the  revised  model)  conformed  to  the  NID  (0»c  )  requirement. 
Additionally,  transformation  of  the  predictor  variables  x  and  -*r  to  yield  a 
model  displaying  no  lack-of-flt,  although  theoretically  possible,  would  ae  a 
long  and  time-consuming  process. 

Fortunately,  the  experimental  design  helped  compensate  for  the  model  lacx- 
of-flt  deficiency.  First  of  all,  the  lack-of-fit  was  comparatively  small. 
Under  optimal  conditions  (lack-of-flt  statistically  Insignificant),  both  the 
mean  square  error  and  the  mean  square  pure  error  estimate  the  population 
variance.  If  there  is  a  lack-of-fit,  the  mean  square  pure  error  estimates  the 
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variance  plus  a  bias  term.  From  the  ANOVA  (regression  and  GLM)  tables,  the 
tabulated  values  for  MSg  and  MSpe  were  54.58  and  52.59  respectively.  Thus, 
although  the  lack-of-flt  was  statistically  significant.  It  was  also  practically 
Insignificant.  In  other  words,  the  development  of  an  alternative  model 
displaying  no  lack-of-flt  would  yield  essentially  (within  1  or  2  51  the  same 
results  as  the  present  model  displaying  lack-of-flt.  Rather  than  Identify  an 
alternate  model  (which  would  not  be  that  much  better  a  predictive  tool  than  the 
model  currently  being  used),  we  decided  to  proceed  with  the  stepwise  model  and 
modify  the  significance  of  the  tests  to  compensate  for  model  lack-of-flt.  Thus, 
all  i's  listed  are  somewhat  higher  than  they  should  be,  and  the  confidence 
Intervals  established  are  somewhat  wider  than  Indicated  In  the  appendix  tables. 

The  final  step  In  determining  model  accuracy  Involved  examining  the 
residual  and  partial  residual  plots  to  verify  the  least-squares  regression 
assumption  that  the  model  errors  were  NID  (0»c^).  Confirming  this  assumption 
confirmed  the  basis  of  the  F  tests  used  to  determine  the  statistical 
significance  of  the  parameters,  and  confirmed  the  statistics  which  produced  the 
tables  listing  confidence  Intervals  of  the  prolactin  concentrations.  A  number 
of  residual  plots  suggested  themselves  Immediately:  residuals  versus  time, 

residuals  versus  predicted  value  of  prolactin  concentration,  and  residuals 
versus  animal  case  number;  studentized  residuals  versus  the  3  plots  just 
mentioned,  and  partial  residual  plots  corrected  for  the  parameters^,  - ^ ,  and 

w  w  • 

4*  ▲ 

Examination  of  the  original  residual  plots  essentially  confirmed  the  NIC 
(0,  -  )  hypothesis.  However,  there  was  one  outlier  in  the  data  set  (case  number 
101,  week  0,  prolactin  concentration  80  ng/ml)  whose  studentlzed  residual  was 
8.16  (Cook's  distance  of  0.097).  This  value  was  most  likely  due  to  an  error  in 
assaying  or  reporting  the  results,  and  was  sufficiently  anomalous  to  be 
discarded  from  the  data  set.  The  residual  plots  were  then  regenerated  and 
rechecked  for  their  distributional  properties.  The  new  plots  (Appendix  F) 
indicated  no  further  problems. 

Diagnostics  to  check  for  model  multlcol 1 inearlty  and  correlation  between 
the  terms  were  then  emoloyec.  Examination  of  the  listed  condition  numbers  and 
matrix  eigenvalues  detected  no  troublesome  values.  This  indicated  that  the 
model  did  not  display  a  significant  degree  of  multlcol 1 inearlty .  Similarly, 
examination  of  the  correlation  matrix  showed  that  correlation  between  the 
estimated  values  of  £  were  all  within  tolerable  limits.  The  highest  degree  of 
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correlation  was  between  the  x  and  the  xz  term*  which  often  occurs  when  using  a 
polynomial  model  In  linear  regression. 

For  future  reference*  and  for  the  sake  of  completeness*  tables  listing 
animal  case  number*  observations  (If  taken)  at  each  week,  predicted  value  of 
prolactin  concentration,  standardized  error  of  prediction,  95%  confidence 
Intervals  on  the  mean  value  of  the  prolactin  concentration,  and  residuals  were 
prepared  as  were  tables  containing  animal  case  number,  regular  and  studentized 
residual  values*  a  graphical  display  of  student  residual  values,  and  influence 
statistics  (such  as  Cook's  0).  These  tables  were  used  to  detect  both  outliers 
and  Influential  data  points  In  the  prolactin  data  set. 

Since  the  null  hypothesis  In  the  study  was  not  rejected  In  the  analysis*  It 
was  necessary  to  determine  the  smallest  difference  between  the  exposed  and  sham- 
exposed  means  that  the  protocol  could  reliably  detect.  A  conservative  estimate 
of  this  sensitivity  was  obtained  by  finding  this  difference  1  in  a  simpler 
experimental  setting.  Since  the  experimental  hypothesis  being  tested  In  this 
alternative  model  was  more  general  than  the  hypotheses  given  In  the  original 
model,  the  difference  obtained  In  the  calculations  would  be  somewhat  larger  than 
the  difference  that  the  ANOVA  design  was  capable  of  detecting. 

To  begin.  It  was  assumed  that  the  experimental  hypothesis  was  merely  one 
testing  the  equality  of  the  means  between  the  exposed  and  sham-exposed  groups 


H  #  s'1 

O'  4  sham-exposed  Z  *  exposed 
Hi :  u  .  r  ,  f  u  . 

i  sham-exposed  exposed 


(3-7) 


This  type  of  hypothesis  could  be  tested  using  a  2-slded  t-test.  The  equation  to 
determine  the  type  II  error  In  this  test  was  then 


This  equation  assumed  equality  In  the  variances  of  the  exposed  and  sham-exposed 


populations. 


;eneral,  this  assumption  was  acceptable  since  there  was  no 


evidence  that  RFR  affected  the  variance  of  prolactin  parameters  differently  in 
the  exposed  and  sham-exposed  groups.  In  equation  (B-8),  the  square  root  of  the 
MSg  estimated  the  population  standard  deviation. 

The  number  of  replications  per  group,  n,  was  computed  by  the  following 
equation: 


and  the  computation  yielded  n  *  63. 

Then#  the  tabulated  value  of  d  was  read  from  the  2-slded  t-test  operating 
curve  for  0.05#  3*  0.10#  and  n=63.  Returning  to  the  original  equation: 

0.45  «  3/7.387  (B-10) 

5  -  3.3242 

Therefore#  the  protocol  was  able  to  detect  a  +  3.32  ng/mL  change  In  resting 
prolactin  concentration  approximately  90%  of  the  time. 

At  the  conclusion  of  the  statistical  analysis#  It  was  evident  that.  If 
there  were  any  RFR-Induced  effects  on  plasma  prolactin  concentration#  these 
effects  were  within  a  range  of  +  3.32  ng/ml  from  the  normal  resting  value. 
Since  this  range  was  within  the  normal  range  of  plasma  prolactin  concentration 
variation  (10  to  15  ng/mL)#  from  a  practical  standpoint#  there  was  no  Indication 
of  RFR-Induced  stress  affecting  animal  resting  plasma  prolactin  concentrations. 

We  gratefully  acknowledge  the  assistance  of  Dr.  Russell  G.  Heikes  of 
Georgia  Tech's  Department  of  Industrial  and  Systems  Engineering  In  developing 
the  statistical  methodology  of  this  appendix. 
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9  CMS  FILEDEF  26  DISK  PROLAC6  LISTING  Al  (BLKSIZE  141  RECFM  VBA  LRECL  133; 

10  CMS  FILEDEF  27  DISK  PROLAC7  LISTING  Al  (BLKSIZE  141  RECFM  VBA  LRECL  133; 

11  ARRAY  WEEK  {33}  WKN3  WKN2  MISSN1  WK0-WK24  MISS25  WKP2  MISS27  MISS28  WKP5 ; 

12  KEEP  X  XSQR  Y  Z  XZ  XSQRZ  CASE; 

13  INFILE  X; 

14  INPUT  CASE  1-3 

15  WKN3  5-6 

16  WKN2  8-9 

17  WKO  11-12 

18  WK1  14-15 

19  WK2  17-18 

20  WK3  20-21 

21  WK4  23-24 

22  WK5  26-27 

23  WK6  29-30 

24  WK7  32-33 

25  WK8  35-36 

26  WK9  38-39 

27  WK10  41-42 

28  WK1 1  44-45 

29  WK12  47-48 

30  WK13  50-51 

31  WK14  53-54 

32  WK15  56-57 

33  WK16  59-60 

34  WK17  62-63 

35  WK18  65-66 

36  WK19  68-69 

37  WK20  71-72 

38  WK21  74-75 

39  WK22  77-78 

40  WK23  80-81 

41  WK24  83-84 

42  WKP2  86-87 

43  WKP5  89-90 

44  ; 

45  MISSN1-. ; 

46  MISS25-.; 

47  MISS27-.  ; 

48  MISS28-. ; 

49  IF  CASE  <  100  THEN  Z  -  0; 

50  IF  CASE  >-  100  THEN  Z  -  1: 
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51  IF  Z  -  1  THEN  CASE  -  CASE  -  100; 

52  DO  I  -  1  TO  33; 

53  X  -  1-4;  XSQR  -  X*X  ;  XZ  -  X*Z;  XSQRZ  -  X*X*Z;  Y  -  WEEK  Cl) ; OUTPUT; 

54  END; 


NOTE:  INFILE  X  IS  FILE  PROLAC  DAT  Al 
NOTE:  126  LINES  WERE  READ  FROM  INFILE  X. 

NOTE:  DATA  SET  WORK.TESTP  HAS  4158  OBSERVATIONS  AND  7  VARIABLES. 

NOTE:  THE  DATA  STATEMENT  USED  0.66  SECONDS  AND  200K. 

55  PROC  CONTENTS; 

NOTE:  THE  PROCEDURE  CONTENTS  USED  0.18  SECONDS  .AND  456K  .AND  PRINTED  PACES  1  TO 

56  PROC  PRINTTO  NEW  UNIT-20; 

NOTE:  THE  PROCEDURE  PRINTTO  USED  0.02  SECONDS  AND  328K. 

57  PROC  SORT  OUT-SCTR; 

58  BY  Z  X  Y; 

NOTE:  DATA  SET  WORK.SCTR  HAS  4158  OBSERVATIONS  AND  7  VARIABLES. 

NOTE:  THE  PROCEDURE  SORT  USED  0.93  SECONDS  AND  295 2K. 

59  PROC  SUMMARY; 

60  BY  Z  X; 

61  VAR  Y; 

62  OUTPUT  OUT-O VLMN  MEAN-MEAN; 

NOTE:  THE  DATA  SET  WORK.OVLMN  HAS  66  OBSERVATIONS  AND  5  VARIABLES. 

NOTE:  THE  PROCEDURE  SUMMARY  USED  0.68  SECONDS  AND  456K. 

63  DATA  S PRO LAC: 

64  SET  SCTR  OVLMN ; 

65  BY  Z; 

NOTE:  DATA  SET  WORK.SPROLAC  HAS  4224  OBSERVATIONS  AND  10  VARIABLES. 

NOTE:  THE  DATA  STATEMENT  USED  0.69  SECONDS  .AND  328K. 

66  PROC  PLOT  NOLEGEND  DATA-SPROLAC ; 

67  BY  Z; 

68  PLOT  MEAN*X- ' X '  Y*X- ’ . ’  /  VAXIS-0  TO  55  BY  5  OVERLAY; 

69  TITLE  'PROLACTIN  SCATTER  DIAGRAM’; 

NOTE:  THE  PROCEDURE  PLOT  USED  1.34  SECONDS  AND  456K  AND  PRINTED  PAGES  3  TO  4. 

70  PROC  PRINTTO  NEW  UNIT-21; 

NOTE:  THE  PROCEDURE  PRINTTO  USED  0.02  SECONDS  .AND  328K. 

71  PROC  PLOT  NOLEGEND  DATA-SPROLAC; 

72  PLOT  MEAN*X- ' X ’  /  VAXIS-0  TO  55  BY  5; 

73  TITLE  'Mean  Plasma  Prolactin  Concentrations  Versus  Time'; 

NOTE:  THE  PROCEDURE  PLOT  USED  1.04  SECONDS  AND  456K  AND  PRINTED  PAGE  5. 

74  PROC  PRINTTO  NEW  UNIT-22; 

75  TITLE  'PROLACTIN  ANALYSIS' ; 

NOTE:  THE  PROCEDURE  PRINTTO  USED  0.02  SECONDS  AND  328K. 
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76  PROC  DATASETS; 


77 

LIST  OF 

MEMBERS 

BEFORE 

UPDATE  OF 

DIRECTORY. 

NAME 

MEMTYPE 

OBS 

TRACKS  PROT 

OVLMN 

/DATA 

66 

1 

SCTR 

/DATA 

4158 

1 

SPROLAC 

/DATA 

4224 

1 

TESTP 

/DATA 

4158 

1 

77 

DELETE 

SCTR; 

78 

DELETE 

OVLMN; 

LIST  OF 

MEMBERS 

AFTER 

UPDATE  OF 

DIRECTORY. 

NAME 

MEMTYPE 

OBS 

TRACKS  PROT 

SPROLAC 

/DATA 

4224 

1 

TESTP 

/DATA 

4158 

1 

NOTE:  THE  PROCEDURE  DATASETS  USED  0.11  SECONDS  AND  456K. 

79  PROC  STEPWISE; 

.80  MODEL  Y  -  X  XSQR  2  XZ  XSQRZ 

81  /  STEPWISE  MAXR; 

NOTE:  THE  PROCEDURE  STEPWISE  USED  0.69  SECONDS  AND  456K  AND  PRINTED  PAGES  6  TO  8. 

82  PROC  PRINTTO  NEW  UNIT-23; 

NOTE:  THE  PROCEDURE  PRINTTO  USED  0.02  SECONDS  AND  328K. 

83  PROC  REG; 

84  MODEL  Y  -  X  XSQR  /  PARTIAL; 

85  ID  CASE; 

NOTE:  ACOV  AND  SPEC  OPTION  ONLY  VALID  WITH  RAWDATA 

NOTE:  THE  PROCEDURE  REG  USED  1.88  SECONDS  AND  648K  AND  PRINTED  PAGES  9  TO  12. 

36  PROC  PRINTTO  NEW  UNIT-24; 

NOTE:  THE  PROCEDURE  PRINTTO  USED  0.02  SECONDS  .AND  328K. 

87  PROC  GLM; 

88  CLASS  X  Z; 

89  MODEL  Y  -  X  Z  X*Z; 

NOTE:  THE  PROCEDURE  GLM  USED  4.13  SECONDS  AND  1032K  AND  PRINTED  PAGES  13  TO  14. 

90  PROC  PRINTTO  NEW  UNIT-25; 

NOTE:  THE  PROCEDURE  PRINTTO  USED  0.03  SECONDS  AND  328K. 

91  PROC  REG; 

92  * - * 

93  * 

94  "  to  obtain  tables  listing  the  variance  inflation  factors.  * 

95  *  influence  statistics,  and  tolerances,  the  following  SAS  •' 


96  "  statements  were  used  in  this  partition:  * 

97  * 

98  *  PROC  REG;  * 

99  *  MODEL  Y  -  X  XSQR  /  TOL  VIF  INFLUENCE;  * 

100  *  ID  CASE;  * 

101  *  OUTPUT  OUT-RPROLAC  P-PREDICT  R-RESID  STUDENT-STUDENT;  * 

102  * 

103  * - *; 


104  MODEL  Y  -  X  XSQR  /  I  SSI  SS2  STB  COVB  CORRB  SEQB  COLLIN 
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105  .  COLLINOINT  ACOV  p  R  CLM; 

106  ID  CASE; 

107  OUTPUT  OUT-RPROLAC  P-PREDICT  R-RESID  STUDENT- STUDENT; 

NOTE:  THE  DATA  SET  WORK.RPROLAC  HAS  4224  OBSERVATIONS  AND  13  VARIABLES. 

NOTE:  THE  PROCEDURE  REG  USED  8.88  SECONDS  AND  648K  AND  PRINTED  PAGES  15  TO  104. 

108  PROC  PRINTTO  NEW  UNIT-26; 

NOTE:  THE  PROCEDURE  PRINTTO  USED  0.02  SECONDS  AND  328K. 

109  PROC  PLOT  DATA-RPROLAC; 

110  PLOT  RES ID*X-’*' ; 

111  PLOT  RES ID’'fPR EDICT- ; 

112  PLOT  STUDENT' *X- ’ * ’ ; 

113  PLOT  STUDENT*PREDICT- ’ * ' ; 

114  TITLE  'PROLACTIN  ANALYSIS' ; 

NOTE:  THE  PROCEDURE  PLOT  USED  1.76  SECONDS  AND  456K  AND  PRINTED  PAGES  105  TO  108 

115  PROC  PRINTTO  NEW  UNIT-27; 

NOTE:  THE  PROCEDURE  PRINTTO  USED  0.02  SECONDS  AND  328K. 

116  PROC  PLOT  DATA-RPROLAC; 

117  BY  2; 

118  PLOT  RESID*CASE» ' * '  /  HAXIS-1  TO  63  BY  2; 

119  PLOT  STUDENT*CASE» ' * '  /  HAXIS-1  TO  63  BY  2; 

120  TITLE  'PROLACTIN  ANALYSIS' ; 

NOTE:  THE  PROCEDURE  PLOT  USED  1.44  SECONDS  AND  456K  AND  PRINTED  PAGES  109  TO  112 

NOTE:  SAS  INSTITUTE  INC. 

SAS  CIRCLE 

PO  BOX  8000 

CaRY ,  N.C.  27511-8000 
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P90LACTIN  ANALYSIS 
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NOTC:  1919  09 S  HAD  MISSING  VALUCS  09  V*9f  OUT  Of  9AWGC  74  MS  HlOOtN 


Residua  Is  versus  animal  ID  number 
(sham-exposure  group). 
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Residuals  versus  animal  ID  number 
(exposure  group). 
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Studentlzed  residuals  versus  time. 
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NOTf:  3272  OSS  HAO  MISSING  VALUCS  691  OSS  HIOOCN 


Studentlzed  residuals  versus  predicted  value  of 
plasma  prolactin  concentration. 
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